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Abstract. One usually expects localized solitons in an integrable system to interact 
trivially. There is an integrable (2+1) -dimensional chiral equation which admits multi- 
soliton solutions with trivial dynamics. This paper describes how to generate explicit 
solutions representing nontrivial soliton interactions: in particular, a head-on collision of 
two solitons resulting in 90° scattering. 
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1. Introduction. 

This paper is concerned with localized solitons in the plane, i.e. in three-dimensional space- 
time. Systems admitting such solitons may be grouped into various distinct classes. For 
example, there are systems admitting topological solitons, the stability of which depends 
on nontrivial topology: these include vortices in the abelian Higgs model [1], and lump 
solutions of sigma-models (with various possible modifications) [2]. These are not solitons 
in the strictest sense; for example, the collision of two solitons is not elastic (some radiation 
is emitted). The feature of topological systems which is relevant here, is that a head-on 
collision results in 90° scattering. In other words, if two solitons approach each other 
along the x-axis and collide at the origin, then two solitons emerge, travelling in opposite 
directions along the y-axis (with slightly less speed, because of the inelasticity). When 
they overlap at the origin, they form a ring rather than a single lump. This phenomenon 
has been observed in numerical simulations [3-7] and understood analytically [8-13]. 

Some other systems that admit localized solitons are those that are completely inte- 
grable. One long-standing example is the Kadomtsev-Petviashvili equation. In this case, 
there are solutions representing n solitons, the interaction of which is completely trivial 
[14]. In particular, if two solitons collide head-on, then they emerge in the same direction 
(without even a phase shift) , and with the same speeds. Another system in which exactly 
the same thing happens is an integrable chiral model [15], which is the subject of this 
paper, and which will be described below (see eqn 1). This system is closely related to one 
introduced in [16], where once again (to quote from that paper) "two-dimensional solitons 
do not interact". One might be tempted to conclude that for integrable systems, local- 
ized solitons on the plane have trivial dynamics. (There are exceptions such as dromion 
solutions of the Davey-Stewartson equations, but these are driven by nontrivial boundary 
conditions at spatial infinity.) 

This picture was, however, undermined by some numerical solutions of the integrable 
chiral equation [17, 7]. These revealed that solitons can also scatter at right angles, as in 
the topological models referred to previously. Solitons in this system have internal degrees 
of freedom, and the numerical results seemed to indicate that solitons can interact either 
trivially or nontrivially, depending on the orientation of these internal parameters. Since 
the system is integrable, one might expect there to be explicit solutions which exhibit 
nontrivial scattering. This paper provides an example of such a solution, and gives a 
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systematic way of constructing many more. 

We must begin, however, by setting up the system. It is defined on (2+1) -dimensional 
space-time, with coordinates (t, x, y) where t denotes time. The chiral field J(t, x, y) is a 
2x2 unitary matrix with det J = 1, and its equation of motion is 



Here d t = d/dt, J t = dJ/dt, etc. The standard SU(2) chiral equation has only the first 
three terms of (1), and not the commutator term: but the latter equation is not integrable. 
By contrast, eqn (1) is a dimensional reduction of the self-dual Yang-Mills equations in 
2+2 dimensions; it is associated with a linear system (see eqn 6), has an inverse scattering 
description [18], passes the Painleve test for integrability [19], and so forth. 

The energy of J is defined to be E = ff £ dxdy, where the energy density £ is 



This quantity E is a positive-definite functional of J, and is conserved [15]. In order to 
ensure finite energy, we impose on J the boundary condition 



where x + iy = rexp(z6>). The matrix J is constant, and J\ is allowed to depend only on 9 
(not on t). 

The simplest nontrivial solution of this system is the static 1-soliton located at the 
origin, which is given by 



Here / is the identity 2x2 matrix, p is the 2-dimensional row vector (1,-2) where z = x + iy, 
p* denotes the complex conjugate transpose of p, and ||p|| 2 = p ■ p* is the norm-squared of 
p. There is numerical evidence that this solution is stable [20]. The energy density goes 
like 8 = 0(r -4 ) as r — > oo, which is the case for all the solutions described in this paper. 



dtiJ^Jt) - d^J^J^ - dyiJ^Jy) + [J~ 1 jy, J ^ J t \ = 0. 



£ = -tr[(J- l J t ) 2 + (J" 1 J,) 2 + (J~ 1 J y ) 2 ] . 



J = J + J 1 (9)r~ 1 + 0(r~ 2 ) asr^oo, 



(2) 
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2. A Scattering Solution. 

In this section, we shall see an example of a solution of (1) representing two solitons which 
undergo 90° scattering. The solution has the form of a product 

'-('-^tX'-tS 2 )- 

where p and q are the row vectors 

V = (M), 

(4) 

q = (l + r 2 ) (l,z)-2ig(z,l), 

with g = t + z 2 . This field J takes values in SU(2); is smooth everywhere on space-time 
(note that p and q are nowhere-zero, which is necessary for smoothness); satisfies the 
boundary condition (2); and satisfies the equation of motion (1) (this follows from the 
construction in the next section). 

The crucial features of this time-dependent solution may be inferred as follows. If 
r 2 = zz is large, then J is close to its asymptotic value Jo, as long as gj ' z 2 = 1+t/z 2 ~ 1. 
But as z approaches ±V~~ t, i.e. as g — > 0, J departs from its asymptotic value: this is 
approximately where the solitons are located. If, therefore, t is negative, then there are 
two solitons, on the x-axis at x ± v / — t; while if t is positive, then there are two solitons 
on the y-axis, at y ss ±y/i. So the picture is of two solitons accelerating towards each 
other, scattering at right angles, and then decelerating as they separate. They accelerate 
as if there were a mutual attractive force, proportional to the inverse cube of the distance 
between them. 

This picture can be confirmed by looking in more detail at the energy density of the 
solution, which is 

s _ 30 1 + 10r2 + 5r4 + 4t2 ( 1 + 2r2 ) - 8t ( x2 - y 2 ) (5) 

[1 + 2r 2 + 5r 4 + At 2 + 8t(x 2 - y 2 )] 2 

Note, first, the symmetry of S under the interchange t h —t, x ^ y, y ^ x: the collision 
is time-symmetric and elastic in this sense, and produces no radiation. For large (positive) 
£, S is peaked, as expected, at two points on the y-axis, namely y ±0.92^. The 
corresponding solitons are not, however, of constant size: their height (the maximum 
value of £) is proportional to l/£, while their radii are proportional to \ft. (The total 



"volume", i.e. the spatial integral of £, is of course constant in time.) In other words, the 
solitons spread out as they move apart. 

Figure 1 illustrates what happens near t = 0. The solitons coming in along the y-axis 
merge to form a ring, and then emerge along the x-axis. 

3. Construction of Solutions. 

This section will indicate how the solution (3) was constructed, and how other solutions 
may be obtained. The procedure is a variant of that in [15]; the latter method was in turn 
pioneered by Zakharov et al [21, 22]. 

The nonlinear equation (1) is the consistency condition for a pair of linear equations 
for a 2 x 2 matrix ip((), where £ is a complex variable (ip also depends on the space-time 
variables t, x, y). This linear pair is 

:= ((d x -d y - W = 

(6) 

Drt := ((dt-Cd y -d x )ij = A 2 ^, 

where A\ and A 2 are 2x2 matrices which are independent of £. Indeed, the integrability 
condition for (6) implies that there exists a J such that 

A^J-^Jt + Jy), A 2 = J~ 1 J x] (7) 

and that this J satisfies the equation of motion (1). Comparing (6) and (7), we see that J 
can be identified with ip(0) _1 . 

The unitarity condition on J follows from an analogous condition on ip(C), namely 

mcmq* = i- (8) 

So the idea is that if we can find a V'(C) such that the unitarity condition (8) holds, and 
such that {Di%lj)ip~ l is independent of ( for i = 1, 2, then J = -(/'(O) -1 is a unitary solution 
of(l). 

The standard way of constructing multi-soliton solutions is to assume that ip has 
simple poles in C, in fact that ip has the form 



where the M k are matrices independent of (. This leads to an n-soliton solution, in which 
the velocity of the kth soliton is determined by the constant complex number fik', one 
consequence is that there is no scattering [15, 16]. 

The scattering solution (3), by contrast, arises from a ip with a double pole (and no 
other poles) in £. In other words, ip is taken to have the form 

^' + cVk^F (10) 

The first requirement on this if) is the unitarity condition (8) . By examining the coefficients 
of the (apparent) poles of •i/K 0^(0*7 one ma y easily show that ip satisfies (8) if and only 
if it factorizes as 

^ (0=|/ "(c-rt n# r-(c-M> 1 



for some 2-vectors p and q. In fact, the same is true if has a pole of order n: the 
unitarity condition is satisfied if and only if ip factorizes into n simple factors of the type 
appearing in (11) (proof by induction on n). 

Next, we have to ensure that (D^ip -1 is independent of £: this imposes differential 
equations on p and q. These are coupled nonlinear equations, and it seems difficult to 
find their general solution. One way of proceeding is to take a solution for the simple-pole 
case (9) with n = 2, of which many are known, and to let ilk — > £t for k = 1,2. If we 
arrange things carefully, this limit gives a solution of the double-pole type (11). (It may 
seem strange that one can take the limit of a family of two-soliton solutions with trivial 
scattering, and obtain a two-soliton solution with highly non-trivial scattering; but in fact 
this is exactly what happens.) 

For the n-soliton solution corresponding to (9), each matrix M k is given in terms of 
a rational meromorphic function fk of one complex variable [15]. (Roughly speaking, fk 
describes the shape of the kth soliton.) In our case, with n = 2, we put ix\ =// + £, 
fj>2 = H — £ and take the limit e — > 0. In order for the resulting ift to be smooth for all 
(t, x, y), it is necessary that f2 — fi — » as e — > 0. So let us write fi = f + sh, f2 = f — sh, 
where / and h are rational functions of one variable. On taking the limit e — > 0, we then 
obtain a ip of the form (11), smooth for all (t,x,y), and satisfying the requirement that 
{Diip)il)~ l be independent of (. Consequently, J = V'(O) -1 is a smooth unitary solution 
of(l). 
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So we have a solution depending on the parameter \i and on the two arbitrary ra- 
tional functions / and h. To get the solution of the previous section, take \i = i (which 
corresponds to the "centre-of-mass" of the system being stationary); this yields (3), with 
p and q being given in terms of f(z) and h(z) by 

P=(l,/), 

(12) 

q = (l + \f\ 2 ) (l,f)-2ig(f,-l), 

where g = tf'(z) + h{z). To get the particular form (4) used in the previous section, one 
takes 

f(z) = z, h(z) = z 2 . (13) 

So the scattering solution belongs to a large family: one may take / and h to be 
any rational meromorphic functions of z. It is not difficult to check that J satisfies its 
boundary condition as r — > oo, irrespective of the choice of / and h. 

The equation (1) does not have rotational symmetry in the xy-plane. So one might 
suspect that the solution (3,4) is somehow exceptional, involving as it does solitons moving 
along the coordinate axes. But in fact one can rotate the whole picture of figure 1 through 
any angle in the xy-plane, by making the choice 

f(z) = z, h(z) = e 2i< V. (14) 

The energy density of the solution corresponding to (14) is the same as in (5), except that 
x 2 — y 2 — $l[z 2 ] is replaced by $l[ex.p(2i(f))z 2 ]. Even though the equation is not rotationally 
symmetric, one can find a new solution which is a rotated version of the original solution. A 
possible way of understanding this is to note that (1) can be written in a Lorentz-invariant 
form, as a Yang-Mills-Higgs equation [23, 20]; so it has a "hidden" Lorentz symmetry. 

4. Concluding Remarks. 

The solution described above, with V'(C) having a double pole, is an example of a 2-uniton. 
The idea of n-unitons was introduced in connection with finding SU(iV) chiral fields on R 2 
[24], and it extends naturally to the corresponding system on R 2+1 [23]. For SU(2) chiral 
fields on R 2 , i.e. the static version of the system of this paper, one uniton is enough (the 
static soliton is a 1-uniton). The original reason for introducing n-unitons was that they 



are needed for static SU(n + 1) chiral fields [24]. What we have seen is that higher unitons 
are also needed for the time-dependent SU(2) case. 

It seems likely that there are many more interesting solutions still to be revealed, both 
in the 2-uniton class discussed in the previous section, and corresponding to higher unitons 
(ip having a higher-order pole in Q. One could, for example, ask what happens when the 
impact parameter of the collision is nonzero; and whether the sizes of interacting solitons 
must necessarily be nonconstant. In addition, it would be interesting to elucidate the role 
of higher unitons in the inverse-scattering [16,18] and algebraic-geometry [23] approaches 
to this system. 
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